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We explore the possibility of using the inhomogeneous magnetic field carried by an Abrikosov 
vortex in a type-II superconductor to localize spin-polarized textures in a nearby magnetic semi- 
conductor quantum well. We show how Zeeman-induced localization induced by a single vortex 
is indeed possible, and use these results to investigate the effect of a periodic vortex array on the 
transport properties of the magnetic semiconductor. In particular, we find an unconventional In- 
teger Quantum Hall regime, and predict directly testable experimental consequences due to the 
presence of the periodic spin polarized structure induced by the superconducting vortex lattice in 
the magnetic semiconductor. 



I. INTRODUCTION 

The possibility of exploring both charge and spin de- 
grees of freedom in carriers has attracted the attention 
of the condensed matter community in recent years. The 
discovery of fcrromagnetism in diluted magnetic semicon- 
ductors (DMS) like GaMnAsi, opened the opportunity 
to explore precisely such independent spin and charge 
manipulations. A fundamental property of diluted mag- 
netic semiconductors (DMS) (both new III-Mn-V and the 
more established II-Mn-VI systems) is that a relatively 
small external magnetic field can cause enormous Zee- 
man splittings of the electronic energy levels, even when 
the material is in the paramagnetic stated. This feature 
can be used in spintronics applications as it allows sepa- 
rating states with different spin. For instance, Fiederling 
et al. had successfully used a II-Mn-VI DMS under effect 
of low magnetic fields in spin- injection experiments—. 

Another utilization of this feature has been discussed 
by Due to the giant Zeeman splitting, a mag- 

netic field with considerable spatial variation can be a 
very effective confining agent for spin polarized carriers 
in these DMS systems. 

Producing a non-uniform magnetic field with nanoscalc 
spatial variation inside DMS systems can be an experi- 
mental challenge. One option is the use of nanomag- 
nets deposited on the top of a DMS layer. In fact, 
nanomagnets have already been used as a source of non- 
homogeneous magnetic field 9 ' 10 . Freire et ali^, for exam- 
ple, have analyzed the case of a normal semiconductor 
in the vicinity of nanomagnets. In such case, the non- 
homogeneous magnetic field modifies the exciton kinetic- 
energy operator and can weakly confine excitons in the 
semiconductor. In contrast, we have examined a DMS in 



the vicinity of nanomagnets with a variety of shape a 5 i 6 i . 
In these cases, the confinement is due the Zeeman in- 
teraction which is hundreds of times stronger than the 
variation of kinetic-energy in DMS. 

Another possibility for obtaining the inhomogeneous 
magnetic fields is the use of superconductors. Nanoscale 
field singularities appear naturally in the vortex phase 
of superconducting (SC) films. Above the lower critical 
field B Cl , in the Abrikosov vortex phase, superconduct- 
ing vortices populate the bulk of the sample, forming 
a flux lattice, each vortex carrying a quantum of mag- 
netic flux 0o/2 = h/(2e). The field of a single vortex is 
non-uniformly distributed around a core of radius r ~ £ 
(where £ is the coherence length) decaying away from its 
maximum value at the vortex center over a length scale 
A (where A is the penetration depth). 

A regular two dimensional electron gas (2DEG) in the 
vicinity of superconductors in a vortex phase has already 
been the subject of both theoretical 1 ^ 12 ! 1 ^ 14 ! 1 ^ 16 ! 1 ^ 18 
and experimental studie»i2*22i2i. In this context, because 
the Lande g factor (and thus the Zeeman interaction) is 
very small in normal semiconductors, the main conse- 
quence of the inhomogeneous magnetic fields is a change 
in the kinetic- energy term of the Hamiltonian. The flux 
tubes do not produce bound states, but act basically as 
scattering center a 14 ! 19 . The Zeeman interaction was ei- 
ther neglected or treated as a small perturbation whose 
main consequence was to broaden the already known Hof- 
stadter butterfly subbands. 

In the work presented here, we consider the case where 
the 2DEG is confined inside a DMS. Due to the giant Zee- 
man interaction in paramagnetic DMSs (whose origin is 
briefly explained in Section II), the Zeeman interaction 
becomes the dominant term in these systems, leading to 
appearance of bound states inside the flux tubes. These 



play a central role in our work and qualitatively change 
the nature of the magnetotransport in these systems, 
compared to the ones previously studied. 

In this article we study a SC film deposited on top of a 
diluted magnetic semiconductor quantum well under the 
influence of an external magnetic field Bq = Bqs, z (see 
Fig. [[J, for all values of Bq between the two possible 
asymptotic limits. For very low applied fields, the SC 
film is populated with few isolated vortices, each vortex 
producing a highly inhomogeneous magnetic field which 
localizes spin polarized states in the DMS. In section ITTT1 
we obtain numerically the energy spectrum and the wave 
functions for the bound states localized by the vortex 
field of an isolated vortex in the DMS QW. 

However, the main advantage of using SC vortices to 
generate confining potentials in the DMSs is the possi- 
bility of varying the distance between them by adjusting 
the external magnetic field. For increasing values of the 
external magnetic field Bq, the vortices in the SC are or- 
ganized in a triangular lattice. The lattice spacing a is 
related to the applied magnetic field by Bq — <j)o / (VSa 2 ) . 
In this limit, the triangular vortex lattice in the SC leads 
to a periodically modulated magnetic field inside the 
neighboring DMS layer. 

Since the magnetic field produced by the SC flux lattice 
creates an effective spin dependent confinement potential 
for the carriers in the DMS, the properties of the super- 
conductors will be reflected in the energy spectrum of 
the DMS. In particular, as the lattice spacing and spa- 
tial dependence of the non-uniform magnetic field (our 
confining potentials) change with Bo, we have a peculiar 
system in the DMS where both the lattice spacing and 
the depth of the potentials are modified by an applied 
magnetic field. 

For low external magnetic fields, the overlap between 
the magnetic fields of independent vortices is small and 
as a consequence, the giant Zeeman effect produces deep 
effective potentials. The trapped states of the isolated 
vortex discussed in section IIIII widen into energy bands 
of spin polarized states. The width of the bands is defined 
by the exponentially small hopping t between neighbor- 
ing trapped states. However, since the flux through each 
unit cell is <p/<j)Q = q/p = 1/2, the energy bands will 
be those of a triangular Hofstadter butterfly^. As long 
as the hopping t is small compared to the spacing be- 
tween consecutive trapped states, this Hofstadter prob- 
lem corresponds to the regime of a dominant periodic 
modulation, which can be treated within a simple tight- 
binding model 22,23 and each band is expected to split into 
p (in this case, p = 2) magnetic subbands. This band- 
structure has unique signatures in the magnetotrans- 
port, as discussed below. Its measurement would pro- 
vide a clear signature of the Hofstadter butterfly, which 
is currently a matter of considerable experimental inter- 

egij 24 i 25 i 26_ 

On the other hand, for a high applied external mag- 
netic field Bq, the magnetic fields of different vortices 
begin to overlap significantly. In this limit, the total 
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FIG. 1: Sketch of the type-II SC - DMS heterostructure in a 
uniform external magnetic field. 



magnetic field in the DMS layer is almost homogeneous. 
Its average is B e z , but it has a small additional periodic 
modulation (the same effect can be achieved at a fixed 
Bq by increasing the distance z between the SC and DMS 
layers). For such quasi- homogeneous magnetic fields, the 
Zeeman interaction can no longer induce trapping; in- 
stead it reverts to its traditional role of lifting the spin 
degeneracy. The small periodic modulation insures the 
fact that the system still corresponds to a 4>/4>q = 1/2 
Hofstadter butterfly, but now in the other asymptotic 
limit, namely that of a weak periodic modulation. In 
this case, each Landau level is expected to split into q 
subbands, with a bandwidth controlled by the amplitude 
of the weak modulation. The case 4>/4>o = 1/2 has q = 1 
and there are no additional gaps in the bandstructure. As 
a result, one expects to see the usual IQHE in magneto- 
transport in this regime. To summarize, as long as there 
is a vortex lattice, the setup corresponds to a <$>/4>q = 1/2 
Hofstadter butterfly irrespective of the value of external 
magnetic field Bq. Instead, Bq controls the amplitude 
of the periodic modulation, from being the large energy 
scale (small Bq) to being a small perturbation (large Bq). 

In section llVl we use a unified theoretical approach to 
analyze how the 2-D modulated magnetic field produced 
by the vortex lattice affects the free carries in the DMS 
QW, and the resulting bandstructures. We obtain the 
energy spectrum going from the asymptotic limit of a 
very small periodic modulation to the asymptotic limit 
of isolated vortices. In the later case, we are able to 
reproduce the results obtained in the Section III. 

As one of the most direct signatures of the Hofstadter 
butterfly, in Section we discuss the fingerprints of the 
band structures obtained in section llVl on the magneto- 
transport properties of these 2DEG, in particular their 
Hall (transversal) conductance. This is shown to change 
significantly as one tunes the external magnetic field be- 
tween the two asymptotic limits. Finally, in Section VI 
we discuss the significance of these results. 
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II. GIANT ZEEMAN EFFECT IN 
PARAMAGNETIC DILUTED MAGNETIC 
SEMICONDUCTORS 

One of the most remarkable properties of the DMS 
is the giant Zeeman effect they exhibit in their param- 
agnetic state, with effective Lande factors of the charge 
carriers on the order of 10 2 - 10 3 . Since this effect plays 
a key role in determining the phenomenology we analyze 
in this work, we briefly review its origin in this section, 
using a simple mean- field picture. 

The exchange interaction of an electron with the im- 
purity spins Si located at positions R^, is: 



interaction of the charge-carriers is, then: 



Hp 



J{r- Ri)S l ■ s. 



(1) 



where J{r) is the exchange interaction. Within a mean- 
field approximation (justified since each carrier interacts 
with many impurity spins) 5j- s — > (Si)s + Si{s), and the 
average exchange energy felt by the electron becomes: 



E„ 



Noxa(S)s 



(2) 



Here, x is the molar fraction of Mn dopants, iVo is 
the number of unit cells per unit volume, (S) is the 
average expectation value of the Mn spins, and a — 
J dru* ^_ o (r)J{r)u c ^^(t 1 ) is an integral over one unit 

cell, with u c g =0 (r) being the periodic part of the con- 
duction band Bloch wave-functions. The usual virtual 
crystal approximation, which averages over all possible 
location of Mn impurities, has been used. Exchange fields 
for holes can be found similarly; they are somewhat dif- 
ferent due to the different valence-band wave-functions. 

In a ferromagnetic DMS, this terms explains the ap- 
pearance of a finite magnetization below Tc- as the Mn 
spins begin to polarize, (S) ^ 0, this exchange interaction 
induces a polarization of the charge carrier spins (s) ^ 0. 
In turn, exchange terms like Si(s) further polarize the 
Mn, until self-consistency is reached. 

In paramagnetic DMS, however, (S) = and charge 
carrier states are spin-degenerate. The spin-degeneracy 
can be lifted if an external magnetic field B is applied. 
Of course, the usual Zeeman interaction —g^ss- B is 
present, although this is typically very small. Much 
more important is the indirect coupling of the charge 
carrier spin to the external magnetic field, mediated by 
the Mn spins. The origin of this is the exchange energy of 
Eq. (J2J, and the fact that in an external magnetic field, 
the impurity spins acquire a finite polarization (S) || B, 
= SBs{gfi B SB/(k B T)), where 5 is the value of the 
impurity spins (5/2 for Mn) and Bs is the Brillouin func- 
tion (for simplicity of notation, we assume the same bare 
(/-factor for both charge carriers and Mn spins; also, here 
we neglect the supplementary contribution coming from 
the Si{s) terms, since typically an impurity spin inter- 
acts with few charge carriers). The total spin-dependent 



7~Lex = ~9^bs- B + N a xa(S)s = g e s^BS- B (3) 



where 



9oS 



-9 



Nqxo: 
fi B B 



SB S 



9HbSB 



(4) 



for electrons, with an equivalent expression for holes. For 
low magnetic fields, g e ff becomes independent of the value 
of B, although it is a function of T and x. Its large 
effective value is primarily due to the strong coupling a 
(large J(r)) between charge carriers and Mn spins. 

In the calculations shown here, we use as DMS pa- 
rameters an effective charge-carrier mass m = 0.5m e and 
9cS — 500, unless otherwise specified. Such values are 
reasonable for holes in several DMS, such as GaMnAs 
and CdMnTe. 



III. ISOLATED VORTEX 

In this section we discuss the case where the free car- 
riers in a narrow DMS quantum well are subjected to 
the magnetic field created by a single SC vortex. This 
situation is relevant for very low applied fields, when the 
density of SC vortices is very low. 

We first need to know the magnetic field induced by 
the SC vortex in the DMS QW. For an isotropic super- 
conductor, this problem was solved by PearlSi. The field 
outside of the SC is the free space solution matching the 
appropriate boundary conditions at the SC surface. Let 
r = \J x 1 + y 2 be the radial distance measured from the 
vortex center, while z is the distance away from the edge 
of the superconductor (here, z is the distance between 
the SC film and the DMS QW, see Fig. [TJ. The radial 
and transversal components of the magnetic field created 
by a single vortex in the DMS layer aro 27 i 28 i 29 i 30 : 



BW{r,z) = 
Bi v) (r,z) = 



»o f°° , „. Ji (kr) exp(~kz - ±£ 2 fc 2 ) 



kdk- 



4ttA 2 J r(k + t) 

J (kr)exp(-kz- ±£ 2 fc 2 ) 



4ttA 2 



kdk- 



r(k + t) 



(5) 
(0) 



4>q = h/e is the quantum of magnetic flux, A and £ are 
the penetration depth and the correlation length of the 
superconductor, r = \/k 2 + A -2 and J U (Q are Bessel 
functions. The term k is a cut-off introduced in 

order to account for the effects of the finite vortex core 
size, which are not included in the London theory 3 ^. 

In Figs. [3 and we show typical magnetic field dis- 
tributions for different ratios of £/A and of z/X. As ex- 
pected, the transversal component is largest under the 
vortex core, and decreases with increasing z (see Fig. 
121- The radial component is zero for r = 0, has a max- 
imum for r ~ £ and then decays fast. This magnetic 
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FIG. 2: The transversal (left) and radial (right) components 
of the magnetic field created by a single vortex in Nb (£ = 35 
nm, A = 40 nm), at different distances z (in units of A). 



Brr 



= <Ao/(4ttA^ 



0.207 T. 



field is qualitatively similar to that created by cylindri- 
cal nanomagnets^. The field distribution depends signif- 
icantly on the properties of the SC. Its maximum value 
is bounded by B max — o /(47rA 2 ). For a fixed A, higher- 
values of B and larger gradients (which are desirable for 
our problem) occur for smaller ratios £/A (see Fig. 
It follows that the ideal SC candidates would be extreme 
type-II (with £/A <C 1), and also have small penetration 
depths A. In order to avoid unnecessary complications 
due to pinning, the superconductor should also have low 
intrinsic pinning (e.g., NbSe^ or MgB2^). Incidentally, 
NbSe2 is also attractive since it was shown that it can 
be deposited via MBE onto GaAs^. In the calculations 
shown here, we use SC parameters characteristic of Nb: 
A = 40 nm, £ = 35nm. 

We now analyze the effects of this magnetic field on 
the DMS charge carriers. Using a parabolic band ap- 
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FIG. 3: The transversal (left) and radial (right) components 
of the magnetic field created by a single vortex at a distance 
z=0.1A, for different values of the coherence length (in units 
of A. Here, A = 40 nm and B max = ^ /(47rA 2 ) ~ 0.207 T. 



proximation, the effective Hamiltonian of a charge car- 
rier inside the DMS QW, in the presence of the magnetic 
field B {v \r, z) of the SC vortex, is [see Eq. ©]: 



2m 



2 -^g cS ^-B^\r,z) (7) 



where m and q are the effective mass and the charge of the 
carrier, and (r, z) is the vector potential, B^ (r, z) = 
V x A^ v '(r, z). For simplicity, we consider a narrow QW, 
so that the motion is effectively two-dimensional. As a 
result, z is just a parameter controlling the value of the 
magnetic field. Generalization to a finite width QW has 
no qualitative effects. 

As discussed in Ref. 0, for a dipole-like magnetic field 
such as the one created by the isolated SC vortex, the 
eigenfunctions have the general structure: 



t/i m (r,<t>) = exp(im<? 



( r ) ex P [*^] 



(8) 



where m is an integer and <fi = tan _1 (?//a;) is the polar 
angle. The radial equations satisfied by the up and down 
spin components can be derived straightforwardly: 



_ld_ 
rdr 

1 d 

rdr 



dr 



HI ~7 / \ 

H — ¥ -gbz(r) - e 



dr 



-I) 2 



-<jb z (r)- 



V'T ( r ) = 9b r (r)ipi(r) 



all lengths are in units of A. The unit of energy is 
E = ft 2 1 (2m\ 2 ) = 0.05 meV if we use m = 0.5m e and 
A = 40 nm. e = E/Eq is the eigenenergy. The magnetic 
fields have been rescaled, B^ v \r,z) = </>o/(47rA 2 ) • b(r) 
(from now on, the distance z between the DMS and 
QW is no longer explicitly specified). Finally, g = 
g cS fi B (bo/{8nX' 2 E ) = g cB /8 if m = 0.5m e . The A^(r) 
terms were left out. This is justified since they are neg- 
ligible compared to the Zeeman term, which is enhanced 
by <7eff ~ 10 2 (this was verified numerically). The bound 
eigenstates e < are found numerically, by expanding 
the up and down spin components in terms of a com- 
plete basis set of functions (cubic B-splines). 

The energies of the ground state (m = 0) and the first 
two excited trapped states (m = ±1) are shown in Fig. 
01 as a function of (a) the ratio £/A, (b) the distance 
z/X to the quantum well, and (c) the ratio g e f[m/m e . As 
expected, the binding energies are largest when the mag- 
netic fields and therefore the Zeeman potential well are 
largest, in the limit z — > and £/A — > 0. Since £/A con- 
trols the spatial extent of the Zeeman trap, the distance 
between the ground and first excited states increase for 
decreasing £/A. All binding energies increase basically 
linearly with increasing g e g. 

The ground-state wavefunction, corresponding to m = 
in Eq. ©, is shown in Fig. Hi. While i/jf\r) ^ be- 
cause of the presence of the radial component, its value 
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FIG. 4: Energy of the ground state (circles) and of the first 
two excited states (triangles), for a charge carrier trapped 
in the Zeeman potential created by an isolated vortex, as a 
function of (a) the coherence length of the SC; here z = 0.1 A, 
m = 0.5m e and g e s = 500; (b) the distance between the 
DMS layer and the SC. Here £/A = 35/40, and the other 
parameters are as in (a); and (c) the value of g c ffm/m e . Here 
£/A = 35/40. In all cases, E = h 2 /(2rn\ 2 ) = 0.05 meV, for 
A = 40 nm. 




FIG. 6: Modulation of the transversal component of the mag- 
netic field Bi(r; z) of the flux lattice in a given area of a type 
II superconductor, for an applied external field Bo = 0.07 T, 
Bo = 0.10 T, B = 0.15 T and B = 0.19 T respectively 
in (a), (b), (c) and (d). As the applied field Bo increases, 
the distance between vortices and the modulation of the total 
field decreases. 

which arc normally incident on the DMS layer, suggest- 
ing possible optical manipulation of these trapped, highly 
spin-polarized charge carriers^. After using quite dras- 
tic analytic simplifications, the results of Ref£ also sug- 
gested that depending on the orientation (sign) of the B z 
component, only one set of excited state (either m > 
or m < 0) is present. Numerically, we find both sets of 
states present, with a small lifting of their degeneracy. 



is significantly smaller than that of the spin-up compo- 
nent -0| ' (r) which is favored by the Zeeman interaction 
with the large transversal magnetic field. The expec- 
tation value of the trapped charge carrier spin has a 
"hedgehog" -like structure, primarily polarized along the 
z-axis, but also having a small radial component (see 
Fig. |3Jd). The general structure of these eigenfunctions 
[Eq. JSJ] has been shown to be responsible for allow- 
ing coupling to only one circular polarization of photons 




FIG. 5: (a)ip\ 0) (r) and of E q- ©. for the ground- 

state, when z = 0.1A(b) The corresponding density of 
charge p(r) — \ipi°\r)\ 2 + \ip^(r)\ 2 , transversal spin den- 
sity s z (r) = 1(1^4° ( r )| 2 — l^i ( r )| 2 ) an£ l radial spin density 
s r (r) =Re(^} 0) >)^ 0) (r)). 



IV. TWO-DIMENSIONAL VORTEX LATTICE 

When placed in a finite external magnetic field B c i < 
Bq < B C 2, the SC creates a finite density of vortices 
arranged in an ordered triangular lattice^. Since each 
unit cell encloses the magnetic flux B a 2 V3/2 = 0o/2 of 
its vortex, the lattice constant a ~ 1/y/Bo is controlled 
by the external field Bq, see Fig. El Consequently, a can 
be varied considerably, depending on the ratio T/Tc of 
the temperature T to the SC critical temperature Tc, 
which sets the value of B c i, and the ratio Bq/B c i. The 
corresponding Hamiltonian for the free carriers in the 
DMS quantum well, in the presence of a SC vortex lattice, 
is given by: 



n 



i 

2m 



p + eA L (f; z) 



1 



g e sVBV ■ B L (r;z) (9) 



Here, — e is the charge of the charge carriers, assumed 
to be electrons. Holes can be treated similarly. We use 
r = (x, y) to describe the 2D position of the charge carrier 
inside the narrow (2D) DMS QW. The magnetic field 
Bl if) created by the triangular vortex lattice is the sum 
of the fields created by single vortices [see Eqs. J5J, JBJ]: 



B L (r;2) = ^BW(f-fl;2) 



R 



B z+^e lG %(z) 



(10) 

where the triangular lattice is defined by R = na(l, 0) + 
mf(l,v3) 5 n, m e Z, and G are the reciprocal lattice 
vectors. The first term is the average field per unit cell, 
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which equals the applied external field Bqz. The second 
term is the periodic field induced by the screening super- 
currents. This term has zero flux through any unit cell 
and decreases rapidly as the distance z between the SC 
and the DMS layers increases. As in the previous sec- 
tion, z is here just a parameter, and we will not write it 
explicitly from now on. 

Likewise, we separate the vector potential in two parts, 
corresponding to each contribution of the magnetic field: 

A L (f) = A (r) + a(r). 

In the Landau gauge, Ao(r) — (0,Bqx,0) and a(r) = 

E^o e<d '%» with a d = li& x Mz)]/\G\ 2 - 

Before constructing the solutions for the full Hamilto- 
nian of Eq. ©, we first briefly review the solutions in 
the presence of only a homogeneous field Bq, in order to 
fix the notation. In this case the Zeeman term lifts the 
spin degeneracy, but it is the orbital coupling that es- 
sentially determines the energy spectrum of the carriers, 
which consists of spin-polarized Landau Levels (LL) : 



(a) 



(b) 



E 



— TlUJr 



1\ 1 

N + -j - -g cS fi B B cr. (11) 



In the Landau gauge the corresponding eigenstates are: 



J N,k v 



r(r) 



-H 



N 



(j + lk y )xa, (12) 



where I — ^Jh/(eB$) is the magnetic length and lo c — 
eBo/m is the cyclotron frequency. N > is the index 
of the LL, iJjv(^) are Hermite polynomials and k y is a 
momentum. x<j are spin eigenstates, a z \a = o~X<n & = 
±1. Each LL is highly degenerate and can accommodate 
up to one electron per 2tt1 2 sample area. The filling factor 
v = n2nl 2 , where n is the 2D electron density in the DMS 
QW, counts how many LL are fully or partially filled. 

This spectrum of highly degenerate LLs is very dif- 
ferent from that arising if only the periodic part of the 
magnetic field and vector potentials are present in the 
Hamiltonian. In this case, the system is invariant to 
the discrete lattice translations. As a result, one finds 
the spectrum to consist of electronic bands defined in a 
Brillouin zone determined by the periodic Zeeman poten- 
tial, the eigenfunctions being regular Bloch states. One 
can roughly think of these states as arising from nearest- 
neighbor hopping between the states trapped under each 
individual vortex, discussed in the previous section. 

In order to construct solutions that include both the 
periodic and the homogeneous part of the magnetic field, 
we need to first consider the symmetries of Hamiltonian 
© • Because of the orbital coupling to the non-periodic 
part of the vector potential Aq, the ordinary lattice trans- 
lation operators T(R) = e^ R ' p do not commute with the 
Hamiltonian. Instead, one needs to define so-called mag- 
netic translation operators^: 



T M (R) 




nib M 



X 



Til a k v 



FIG. 7: (a) Rectangular magnetic unit cell containing 2 vor- 
tices, and its lattice vectors, (b) The associated magnetic 
Brillouin zone and the location of the special, high-symmetry 
points r, X and M. 



It is straightforward to verify that these operators com- 
mute with the Hamiltonian, \H, Tm{R)\ — 0, and that 



T m (R)Tm{R') = e 



-jt B 0R':cRy 



T M {R + R'). 



It follows that these operators form an Abelian group 
provided that we define a magnetic unit cell so that the 
magnetic flux through it is an integer multiple of 4>q. In 
our case, as already discussed, the magnetic flux through 
the unit cell of the vortex lattice is precisely </)q/2. We 
therefore define the magnetic unit cell to be twice the 
size of the original one. We use a rectangular unit cell, 
as shown in Fig.[7^. The new lattice vectors are a = (a, 0) 
and b = (0,6), with b = ai/3, and R nrn — na + mb. The 
basis consists of two vortices, one placed at the origin 
and one placed at S = (a + b)/2. The associated mag- 
netic Brillouin zone is k x € (— 7r/a, tt/o\; k y G (— n/b, ir/b] 

and the reciprocal magnetic lattice vectors are G n ^ m = 
n(2Tt/a)x + m(2n/b)y, n,m £ Z (see Fig. 0d). 

With this choice, the wave-functions must also be 
eigenstates of the magnetic translation operators: 



T M (R)^ k (r) = e^~ k (r). 



(13) 



We thus need to expand the eigenstates in a complete 
basis set of wavefunctions which satisfy Eq. I|13|) . Such a 
basis can be constructed from the LL eigenstates, since 
the large degeneracy of each LL allows one to construct 
linear combinations satisfying Eq. (|13|! 3b : 



1 



E 



Vfc H + ^n, CT ( f ) ( 14 ) 



where Nt is the number of magnetic unit cells and k is 
a wavevector in the magnetic Brillouin zone. 

As a result, we search for eigenstates of Hamiltonian 
@ of the general form: 



(15) 
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Here, d/v CT (fc) are complex coefficients characterizing the 
contribution of states from various LLs to the true eigen- 
states. Spin-mixing is necessary because [H, <f z ] ^ 0. 

Note that since the periodic potential may be large 
(due to the large Lande factor) we cannot make the cus- 
tomary assumption that it is much smaller than the cy- 
clotron frequency, and thus assume that there is no LL 
m bd n g36 i 37 i 38 J 39. Instead, we mix a large number of LLs, 
so that we can find the exact solutions even in the case 
when the cyclotron frequency is much smaller than the 
amplitude of the periodic potential. 

The problem is now reduced to finding the coeffi- 
cients dpjaik), constrained by the normalization condi- 
tion J2n a \dNa{k)\ 2 = 1. The Schrodinger equation re- 
duces to a system of linear equations: 



[E% -EN,a]dN<r(k)- 



g c ff^B 



(16) 



N 



where 



dm* N% {f) B L (r;z)-B z * N , g(r) 



Bq{z) exp[~i(k x G y - k y G x )l 2 }I§ N , 



Here (see Ref. |3 



1 N,N' 



(G) 



G x — iG y 
G 



N — N 



where m = min(A, A'), M = max(7V,7V'), G = l 2 G 2 /2. 
L™(x) are associated Laguerre polynomials. The Fourier 
components Bg(z) of the magnetic field can be calculated 
straightforwardly, see Eqs. JSJ, © and ijTU|) . For the 
reciprocal vectors of the magnetic unit cell, we find that 
(z) — if n + m is an even number, otherwise: 



B G 



BMz) = (-iG x ,-iG y ,G) 



Bne 



-Gz-±G 2 £ 2 



Gt(G\ + r) 



where r = V G 2 A 2 + 1. 

Eq. (1 1 (if) neglects the periodic terms proportional to 
a(f ; z) . These are small compared to the periodic terms 
related to B^ (r ; z) , since the latter are multiplied by 
g c ff 3> 1 (we verified this explicitly). We numerically 
solve Eq. (|16|l . typically mixing LL up to TV = 40 
and truncating the sum over reciprocal lattice vectors in 
6 at n>, to the shortest 1600. These cutoffs are such that 



the lowest bands eigenstates tp^ (?) and dispersion E~ 
(a is the band index) are converged and do not change 
if more LL and/or reciprocal lattice vectors are included 
in the calculation. 

A first test for our numerical results is to compare re- 
sults obtained in the lattice limit Bq — > 0, a — > oo, with 




FIG. 8: The energy spectrum in the lattice case with a — 
8.5A (solid lines) compared to lowest energy eigenstates of the 
isolated vortex (dashed lines). SC parameters correspond to 
Nb and z = 0.1 A. 



the energy spectrum of the isolated vortices, obtained 
in the previous section. One such comparison is shown 
in Fig. 03 where the band dispersion for a lattice with 
a = 8.5A is shown to have the location of the lowest bands 
in good agreement with the eigenenergies of the isolated 
vortex. The agreement improves for larger a values. 

Another check on our results is to investigate the dis- 
persion of the lowest-energy bands, in this limit. As dis- 
cussed before, in the absence of the orbital coupling to 
the Aq component of the vector potential, one expects 
a simple tight-binding dispersion for the lowest bands, 
with an effective hopping t characterizing the overlap 
between eigenfunctions trapped under neighboring vor- 
texes. In the presence of the orbital coupling, the hop- 
ping matrices pick up an additional phase phase factor 
proportional to the enclosed flux (the Peierls prescrip- 
tion). This is responsible for lifting the degeneracy of 
each tight-binding band. The resulting sub-band struc- 
ture, and in particular the number of sub-gaps opened 
in each tight-binding band, are known to depend only on 
the ratio 0/</>o> where 4> is the flux of the applied magnetic 
field through the unit cell of the periodic potential. In 
fact, for 0/0o = p/li where p and q are mutually prime 
integers, each tight-binding band splits into precisely q 
subbands^. This problem has been well studied as one 
of the asymptotic limits of the Hofstadter butterfly, cor- 
responding to a large periodic modulation and a small 
applied magnetic field. 

In our case, 4>/4>q = 1/2 and we expect each tight- 
binding band to split into two sub-bands. This is indeed 
verified in all the asymptotic limits a 3> A where the 
tight-binding approximation is appropriate (in some of 
the plots that we show, such as Fig. |S| this splitting is too 
small for the lowest band and is not visible on this scale) . 
In fact, we can even fit the dispersion of the bands, in this 
asymptotic limit. For a triangular lattice with nearest 
neighbor hopping t (real number), it is straightforward 
to show that the dispersion when a magnetic field with 



— • — • - 

t (meV) 
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FIG. 9: Fit of the energy spectrum of the two lowest energy 
sub-bands (circles) with the expressions of Eq. 1171 appropri- 
ate in the asymptotic limit a>A (lines). This energy spec- 
trum corresponds B x — B y = (see text) and a = 4A, and it 
is shown as a function of k x , with k y — 0. Its overall energy 
values have been shifted for convenience. The inset shows the 
effective hopping t extracted from this fit, which decreases 
exponentially with increasing distance between vortices. 



1/2 is added is (see Ref. EJ: 



E(k x , k y ) — ±2t-J 1 + cos 2 (k x a) — cos(k x a) cos(k y b). 

(17) 

In our case, the hopping t between states trapped un- 
der neighboring vortices is not a real number, even if we 
set A (r) — . The reason is that the wave- functions 
have a non-trivial spinor structure (see Eq. ©) which 
leads to a complex value of t (t is just a matrix el- 
ement related to the overlap of neighboring wavefunc- 
tions) . To avoid complications coming from dealing with 
the phase of t and the changes induced by it on the sim- 
ple dispersion of Eq. (|17(l . we test the fit for a "toy 
model" in which we set the in-plane magnetic field to 
zero: B x (x,y) = B y (x.y) = 0. In this case, the spin is a 
good quantum number, the ground-state wavefunctions 
are simple s-type waves, the corresponding hopping t is 
real and Eq. (|17fl holds. We show a fit for such a case in 
Fig. along the k y — line in the Brillouin zone. This 
fit allows us to extract a value for the effective hopping 
t. As expected (see inset of Fig. |§J), t decreases expo- 
nentially with increasing distance a between neighboring 
vortices. 

The presence of the in-plane components changes the 
structure of the wavefunctions and leads to complex val- 
ues of t, modifying the dispersion from the simple Eq. 
(|17fl form. Indeed, as one can see from Fig. 03 the dis- 
persion is now quite different in shape than the one shown 
for the toy model in Fig. |5] In fact, although there are 
two distinct subbands E(k x , k y ) for any point in the Bril- 
louin zone, their gaps do not overlap and so there is no 
true subgap appearing in the band-structure. 
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FIG. 10: Band structure for the carriers in the DMS QW for 
Bo of (a) 0.07T;(b) 0.10T; (a) 0.15T;(a) 0.19T. The gaps are 
marked by shaded regions, and each is labelled with its Chern 
number (see text). 



These results, corresponding to the asymptotic limit 
of a large modulation and small applied field, show that 
our numerical scheme based on expansion in terms of 
multiple LLs is working well even in this most unfavor- 
able limit. The other asymptotic limit where our results 
can be easily verified against known predictions is the 
limit of a large applied field and small periodic modula- 
tion. In this case, the small modulation is expected to 
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lift the degeneracy of each LL, but not to lead to mixing 
amongst the LL. This problem has also been studied ex- 
tensively22i2£i21 and the resulting spectrum is also known 
to depend only on the ratio 4>/4>o — p/q. Unlike in the 
tight-binding limit, here each LL splits into p subbands. 

In our case, p = 1 and therefore we expect no supple- 
mentary structure in the LLs. This is indeed verified, as 
shown, for example, in Fig. 1101 where we plot the evolu- 
tion of the electronic band structure as a is varied. For 
a 3> A (panel (a)) we see the emergence of the tight- 
binding structure discussed above. As a decreases with 
increasing Bq (panel (d)), we indeed see the emergence of 
nearly equidistant Landau levels, which still exhibit some 
dispersion due to the weak periodic potential. Because 
of the large, quasi-uniform Zeeman interaction in this 
limit, all these states are mostly spin- up and the splitting 
between consecutive bands corresponds to the cyclotron 
energy. The intermediary cases (panels (b) and (c)) cor- 
respond to situations where neither asymptotic limit is 
appropriate. In such cases one needs to perform nu- 
merical simulations to find the resulting band-structure. 
The most direct signature of this strongly field-dependent 
band-structure is obtained in magnetotransport measure- 
ments, which we proceed to discuss now. 



V. INTEGER QUANTUM HALL EFFECT 

Following the discovery of the integer quantum Hall 
effect in a two dimensional electron gas (2DEG) in a 
strong magnetic field, Laughlin demonstrated that the 
Hall conductance of a non interacting 2DEG is a mul- 
tiple of e 2 jh if the Fermi energy lies in a mobility gap 
between two Ll4i. Later, Thouless et al. argued that 
the quantization of the Hall conductivity a xy in period- 
ically modulated systems has a topological nature and 
therefore it occurs whenever the Fermi energy lies in an 
energy gap, even if the gap lies within a Landau levels. 
Using the Thouless formula, which we discuss below, the 
Hall conductances corresponding to various types of lat- 
tices, in either of the two asymptotic limits of strong or 
weak periodic potential, have then been investigated^. 

In our case, it is necessary to use a more general 
method to calculate the Hall conductance. We briefly 
review it here. This approach is inspired by the work of 
Kohmoto^ and Uso\*2S, who showed that it is possible 
to use the topology of the band-structure to calculate, 
in a relatively simple way, the contribution to the Hall 
conductance of each filled subband. 

When the Fermi level is inside a gap, the total Hall 
conductivity equals the sum of the contributions from all 
the fully occupied bands: a xy = J2 a a ^'v > ®(^F ~ )• 
The contribution of each occupied band is2&: 



4$ = * 



e 

2nh 



dk / df 



cm- 

k 



(a): 



dui a) 



dui a) * 



dui a) 



where the integrals are carried over the magnetic Bril- 
louin zone and the magnetic unit cell, respectively. Here, 



tpj?' (r)e 



ik-f 



is the Bloch part of the band 



eigenstate. Using Eqs. (TH|l and lfH)|l we can perform the 
real space integrals explicitly to obtain: 



<7<«> = - 
h 



1 + — Im / dk^ 

27T / ^ 

J N<7 



(fe) 



dk v 



dk x 



The first term is the unit conductance contribution of 
each LL (band) , expected in the absence of the periodic 
modulation - the usual IQHE. The second term can be 
rewritten as 



dke 



Vj: x A {a \k) 



where 



A {a \k) 



Na 



(20) 



(21) 



This term can be calculated using Kohmoto's argu- 
ments 43 . The magnetic Brilloum zone is a T 2 torus and 
has no boundaries. As a result, if the gauge field A^ a \k) 
is uniquely defined everywhere, Stoke's theorem shows 
that this integral vanishes. The gauge field is related to 
the global phase of the band eigenstates: if ip^ (r) — > 

e i/(£)^ Q )(f) 5 then A^{k) -> A^{k) - iVj:f(k). It fol- 
lows that the gauge field can be uniquely defined if the 
phase of any one of the components dS a [k) (and thus, the 
global phase of the band cigenfunction) can be uniquely 
defined in the entire magnetic Brillouin zone. This is pos- 
sible only if the chosen component d^ a {k) has no zeros 
inside the magnetic Brillouin zone, in which case we can 
fix the global phase by requesting that this component be 
real everywhere. In this case, as discussed, a x °^ — e 2 /h. 

If there is at least one point ko where d^.(ko) = 0, in its 
vicinity the global phase must be defined from the condi- 
tion that some other component d^)^, (fc), which is finite 
in this region, is real. The definitions of the global phase 
inside and outside this vicinity of k$ are related through 
a gauge transformation; moreover, the torus now has a 
boundary separating the two areas. Applying Stokes' 
theorem, it follows immediately" 4 ^ that Aeriy' = \Sq. 
The integer So is the winding number in the phase of 
d^) ,(k) (when dj^(fe) is real) on any contour surround- 
ing kg. If d^(k) has several zeros inside the Brillouin 
zone, then one has to sum the winding numbers associ- 
ated with each zero: 



(a) 



(22) 



(18) 



Thus, once the coefficients d^ a {k) are known, one can 
immediately find the contribution of the band a to a xy 
by studying their zeros and their vorticities. 
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FIG. 11: Values in the magnetic Brillouin zone of the d%'„(k) 



coefficients corresponding to the second band a = 2, in the 
band-structure obtained for Bo = 0.095T, z = 8 nm and 
g = 500. a) d^l(k) is chosen to be real everuwhere. It has 
zeros in the corners and center of the magnetic Brillouin zone. 



b) |d3 2 l(fe)| has maxima where d^Uk) has zeroes 



>' 2 h 



As an illustration of this method, we calculate the con- 
tribution of the second band a = 2, for B n = 0.095T, 
z = 8 nm and g = 500. We first choose the overall phase- 

(2) — * 

factor for the corresponding coefficients d y N (k) defining 
this band so that one of them (specifically, here we chose 
N = 2) is real throughout the Brillouin zone, see Fig. 
ITTl We see that this component has zeros in some high- 
symmetry points, signalling a potentially non-trivial con- 
tribution to a xy . In order to find the corresponding vor- 
ticities, we investigate any other component that has no 

(2) 

zeros at the positions of the singular points of rf^j. In 
general, this other component will have complex values. 
In the lower panel of Fig. ^2 we plot the absolute value 
of dg 2 ] , which is indeed finite at all zeros of d^l . 

The vorticities can be found by investigating the vari- 

(2) (2) 

ation in the phase of d 3 j around the zeros of d\ i. The 



2,r 



phase-map of G? 3 2 j is shown in Fig. I|12f) . where red repre- 
sents a phase of tt and blue a phase of — 7r. Any bound- 
ary red-blue indicated a non-zero winding number. These 
winding numbers S m (in units of 2ir) are the phase incur- 

(2) 

sions of <fg j upon going anticlockwise around the singular 
points. In this case, it is clear that the phase incursions 
are of — 2tt. As there are two singular points in the MBZ 
(the corners count as 1/4 each), the conductance of the 

subband 2 is ci 2 ' = — 2) = — (x)- The same over- 

all value is obtained by looking at the vorticities of any 



(2) 

2,r 



other wavefunction that is finite at the zeros of d 

Larger winding numbers are also possible. As a second 
example, in Fig. 1131 we show the phase map for the com- 
ponent dg 6 j for the sixth band. In this case, we requested 

that the component d^ be real, and we found its zeros to 
be again in the center and corners of the Brillouin zone. 
Now, the phase incursion about the zeros are — 2 x 2tt 




FIG. 12: Phase map of dg 2 ' inside the MBZ. Blue represents 
—n and red represents 7r. The phase winds by — 2-zr around 
both the center and the corner points. 



and therefore the conductance of the 6 th subband is, in 

this case, axy = ~ 4) = — ^(x)- These integers 

are called Chern numbers. 

As discussed by Avron et alM>, based on the topologi- 
cal nature of the Hall conductance it can be shown that 
the Chern number of a subband does not change unless 
the subband merges with a neighboring subband. In this 
case, the Chern number of the newly formed band equals 
the sum of the Chern numbers of the two subbands. Simi- 
larly, if a band splits into one or more subbands, its Chern 
number is re-distributed over the two or more subbands. 
In the system discussed here, a variation of the applied 
magnetic field Bq changes the band structure, opening 
and closing energy gaps. By calculating the energy spec- 
trum of the system for a given range of the magnetic 
field, we map the opening and closing of the gaps be- 
tween the few lowest mini-bands as a function of Bq (see 
Fig. HT)jl . We than calculate the Hall conductance of each 
of the lowest seven mini-band for each configuration of 
gaps using this method. 

It is important to point out that we recover the ex- 
pected values of the Hall conductance of the minibands 
in both limits of strong and weak modulation for a tri- 
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FIG. 13: Same as in Fig. I12l but for component d^l(k). 
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FIG. 14: Density of states and hall conductivity as a function 
of the Fermi Energy for two different applied magnetic fields 
B . 



angular lattice and p/q = 1/2. In Fig. ^3 we see in 
panel (b) that the first two jumps in the conductivity are 
le 2 /h and then — le 2 //i, as expected for the p — 2 sub- 
bands of tight-binding limit of the triangular Hofstadter 
butterfly 40 . In the limit of weak modulation, we obtain 
the usual IQHE, as expected, since p = 1. 

Such a calculation allows us to predict the sequence of 
plateaus in the IQHE if one keeps Bq fixed and varies the 
concentration of electrons in the DMS (for instance, by 
varying the value of a back-gate voltage). Such predic- 
tions are shown in Fig. 1141 where we plot the density of 
states (lower panels) and the corresponding Hall conduc- 
tivity (upper panels) as a function of the Fermi energy 
for band-structure corresponding to two different values 
of Bq. Of course, one needs disorder in order to observe 
the IQHE. We did not include disorder in this calcula- 
tion, but we know from Laughlin's arguments that the 
value of the Chern numbers remains the same, in its pres- 
ence. The plots in the upper pannels of Fig. ^]are thus 
only sketches of the expected Hall conductances, in these 
cases. The main observation is that the Hall conductiv- 
ity does not increase monotonically as a function of Ep, 
as is the case in the usual IQHE. This, of course, is due 
to the presence of the periodic potential induced by the 
vortex lattice, and clearly signals the formation of the 
Hofstadter butterfly. 

We can summarize more efficiently the information 
shown in plots like Fig. [21 in the following way. The 
two important parameters are the values of the electron 
concentration when the Fermi level is in a gap, and the 
value of a xy expected for that gap. Such a plot is shown 
in Fig. 1151 The various colored regions are centered on 
the values of the electron concentrations where the Fermi 
energy should be in a gap. These regions are assigned an 
arbitrary width (calculations involving realistic disorder 
are needed to find the widths of these plateaus) and are 
labelled with the integer i defining the quantized value 
of <J xy = i(e 2 /h). As the external magnetic field is var- 




012 B(T) 0.14 



FIG. 15: conductivity a xy in units of e 2 /h, as a function of the 
magnetic field B and the charge carrier density n. The colored 
areas mark the gaps between the bands. For some ranges 
of B, neighboring bands touch and some of the gaps close. 
The integers give the quantized values of the cr xy plateaus. 
For example, er : 
n= 5.5 x 10 10 c 
quantized every time the Fermi level is inside a gap. The 
parameters are the same as for Fig. HJ 



as a function of B, at a constant density 
~ 2 (red line) is shown in the inset. It is 



ied, some subbands merge and then separate, changing 
their individual Chern numbers in the process. This plot 
allows us to predict the sequence of Hall plateaus for a 
constant value of the electron density n and varying ex- 
ternal field, as shown in the inset for n — 5.5 x 10 10 
cm~ 2 . The non-monotonic sequences of plateaus as -Bo 
decreases clearly signals the appearance of tight-binding 
bands, which in turn show that the Zeeman potential 
is strong enough to localize spin-polarized charge carries 
under individual vortices. 



VI. CONCLUDING REMARKS 

In this paper we presented detailed numerical and ana- 
lytical calculations aimed at investigating the novel spin 
and charge properties of a magnetic semiconductor quan- 
tum well in close proximity of a superconducting flux line 
lattice. First we have shown how the single supercon- 
ducting vortex localizes spin polarized states, with bind- 
ing energies within the accessible range of several local 
spectroscopic probes, as well as transport measurements. 
We then turned our attention to the case of a periodic 
flux line lattice, and presented the results of a numeri- 
cal framework able to interpolate between the inhomoge- 
neous (low) field regime of dilute vortex lattice and the 
homogeneous (high) field regime characterized by Lan- 
dau level quantization. Our numerical scheme not only 
reproduces the energy spectrum of the isolated vortex 
limit, but the spin-polarized electronic band structure we 
obtain within this framework also matches with the an- 
alytical tight binding calculations applied for the dilute 
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vortex lattice limit as well. Between the two extreme 
field limits we investigated the momentum-space topol- 
ogy of the Bloch wave functions associated with the spin 
polarized bands. By using the connection between the 
wave function topology and quantum Hall conductance, 
we showed how the consequences of the 1/2 Hofstadtcr 
butterfly spectrum can lead to experimentally observ- 
able effects, such as a non-monotonic 'staircase' of Hall 
plateaus as they appear under a varying magnetic field 
or carrier concentration. We paid special attention to 
provide realistic systems and materials parameters for 
each experimental configuration we suggested. All indi- 
cations from our theory seem to suggest that the mag- 
netic semiconductor-superconductor hybrids can be fab- 
ricated with presently available molecular-beam epitaxy 



will provide us with a rich variety of transport and spec- 
troscopic phenomena. 
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